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APPLICATIONS 
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Abstract. Consider a mean-reverting equation, generalized in the sense it is 
driven by a 1-dimensional centered Gaussian process with Holder continuous 
paths on [0, T] (T > 0). 

Taking that equation in rough paths sense only gives local existence of the 
solution because the non-explosion condition is not satisfied in general. Under 
natural assumptions, by using specific methods, we show the global existence 
and uniqueness of the solution, its integrability, the continuity of the associ- 
ated Ito map and we provide an L p -converging approximation with a rate of 
convergence (p 1). The continuity of the Ito map ensures a large deviation 
principle for that generalized mean-reverting equation. 
Finally, we study a generalized mean-reverting pharmacokinetic model. 
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1. Introduction 

Let W be a 1-dimensional centered Gaussian process with a-H61der continuous 
paths on [0, T] (T > and a e]0, 1]). 

Consider the stochastic differential equation (SDE) : 

(1) X t =x + [ (a- bX u ) du + a f X£dW u ; t G [0, T] 

Jo Jo 

where, x > is a deterministic initial condition, a,b,a > are deterministic 

constants and (3 satisfies the following assumption : 

Key words and phrases. Rough paths, rough differential equations, large deviation principle, 
mean-reversion, Gaussian processes. 
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Assumption 1.1. The exponent (3 satisfies : /3 €]1 — a, 1]. 

When the driving signal is a standard Brownian motion, equation taken in the 
sense of Ito, is used in many applications. For example, it is studied and applied 
in finance by J-P. Fouque et al. in [5] for f3 £ [1/2, 1[. The cornerstone of their 
approach is the Markov property of diffusion processes. In particular, their proof 
of the global existence and uniqueness of the solution at Appendix A involves S. 
Karlin and H.M. Taylor [S], Lemma 6.1(h). Still for f3 £ [1/2, 1[, the convergence of 
the Euler approximation is proved by X. Mao et al. in |15) . For j3 ^ 1, equation ([!]) 
is studied by F. Wu et al. in [19]. Recently, in [16_, N. Tien Dung got an expression 
and shown the Malliavin's differentiability of a class of fractional geometric mean- 
reverting processes. 

Equation (JTJ is a generalization of the mean-reverting equation. In this paper, 
we study various properties of (|T|) by taking it in the sense of rough paths (cf. T. 
Lyons and Z. Qian [14 ). Note that Doss-Sussman's method could also be used 
since Q is a 1-dimensional equation (cf. H. Doss [3] and H.J. Sussman [IB])- A 
priori, even in these senses, equation (|T|) admits only a local solution because it 
doesn't satisfy the non-explosion condition of [7J, Exercice 10.56. 

At Section 2, we state useful results on rough differential equations (RDEs) and 
Gaussian rough paths coming from P. Friz and N. Victoir [7J. Section 3 is devoted 
to study deterministic properties of . We show existence and uniqueness of the 
solution for equation 0, provide an explicit upper-bound for that solution and 
study the continuity of the associated Ito map. We also provide a converging ap- 
proximation with a rate of convergence. Section 4 is devoted to study probabilistic 
properties of ([lj) ; basic properties, various integrability results and a large deviation 
principle. Finally, at Section 5, we study a pharmacokinetic model based on a par- 
ticular generalized mean-reverting (M-R) equation (inspired by K. Kalogeropoulos 
et al. HTJ). 

2. Rough differential equations and Gaussian rough paths 

Essentially inspired by P. Friz and N. Victoir [7J, this section provides useful defi- 
nitions and results on RDEs and Gaussian rough paths. 

In a sake of completeness, results on rough differential equations are stated in 
the multidimensional case. 

In the sequel, |.| denotes the euclidean norm on R d and ||.||a4 the usual norm 
onMj(M) (d£W). 

Consider Dt the set of subdivisions for [0, T] and 

A T = {(s,t) £ R 2 + : sC s < t sC T) . 
Let T N (R d ) be the step- TV tensor algebra over R d (N £ N*) : 

N 

T N (R d ) =($(R d f\ 

i=0 

For i = 0, . . . , N, (R d )® 1 is equipped with its euclidean norm \\.\\i, and the canonical 
projection on (R d )® 1 for any Y £ T N {R d ) is denoted by Y\ 

First, let's remind definitions of p-variation and a-H61der norms (p ^ 1 and a £ 
[0,1]): 
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Definition 2.1. Consider y : [0, T] -> R d : 

(1) T/ie function y has finite p-variation if and only if 

(\D\-1 \ 1/P 
\\y\\p-var;T = SU P ^ WVrk+i ~ Vr JN < °°' 



D={r k }GD T 



fe=l 



(2) The function y is a-Holder continuous if and only if 

ii I, \\yt~y s \\ ^ 

\\y\\a-Hdl;T = SUp — — < 00. 

( s ,()eA T |i- s l 

In the sequel, the space of continuous functions with finite p-variation will be de- 
noted by : 

C p " var ([0,T];R d ) . 
The space of a-Holder continuous functions will be denoted by : 

C Q - H61 ([0,T];R d ) . 

If it is not specified, these spaces will always be equipped with norms ||.|| p _var;T and 
||-||a-Hoi;T respectively. 

Remark. Note that : 

CQ -H61 ([ 0)T ]. R d) c C-l/a-var ([ 0)T ]. R ^ _ 

Definition 2.2. Let y : [0, T] — > R d be a continuous function of finite l-variation. 
The step-N signature of y is the functional Sn(v) '■ At — > T N (M. d ) such that for 
every (s, t) € At and i = 1, . . . , N, 

S%- S ,t{y) = 1 and S l N . Sjt (y) = / dy ri ® • • • <g> dy u . 

J s<r 1 <r 2 <-<r i <t 

Moreover, 

G N (R d ) = {S N ., , T (y);y e C 1 - var ([0,T]:R d )} 
is the step-N free nilpotent group over R d . 

Definition 2.3. A map Y : At — > G N (R d ) is of finite p-variation if and only if 

(\D\-1 \ 1/P 

\\Y\\ p . var . tT = sup V ||y rji , rfc+1 ||£ <°° 

D={r k }eD T y k=1 j 

where, ||.|| c is the Carnot-Caratheodory's norm such that for every g € G N (R d ), 

\\g\\ c = inf {length(y);y e C 1 ' var ([0,T];M d ) and S N;0 , T (y) = g) . 

In the sequel, the space of continuous functions from At into G N (R d ) with finite 
p-variation will be denoted by : 

C p - var ([0,T];G w (]R d )). 

If it is not specified, that space will always be equipped with || .|| p - V ar;T- 

Let's define the Lipschitz regularity in the sense of Stein : 

Definition 2.4. Consider 7 > 0. A map V : R d — > R is ^-Lipschitz (in the sense 
of Stein) if and only if V is C^J on R d , bounded, with bounded derivatives and 
such that the \j\-th derivative of V is {^-Holder continuous is the largest 
integer strictly smaller than 7 and {7} = 7 — 

The least bound is denoted by \\V\\u p -i. The map \\.\\ii p i is a norm on the vector 
space of collections of 7- Lipschitz vector fields on R d , denoted by Lip 11 
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In the sequel, Lip 7 (E d ) will always be equipped with ||.||ii p -y. 

Let w : [0, T] — > M. d be a continuous function of finite p- variation such that a 
geometric p-rough path W exists over it. In other words, there exists an approxi- 
mating sequence (w n , n £ N) of functions of finite 1-variation such that : 

lim d p . var;T [S [p] 0");W] =0. 



We remind that if V — (Vi, . . . , Vd) is a collection of Lipschitz continuous vector 
fields on R d , the ordinary differential equation dy = V(y)dw n , with initial condition 
I/O € R i admits a unique solution. 

That solution is denoted by irv(0,yo;w n ). 

Rigorously, a RDE's solution is defined as follow (cf. [7], Definition 10.17) : 

Definition 2.5. A continuous function y : [0,T] — > M. d is a solution of dy — 
V(y)dW with initial condition yo £ M. d if and only if, 

lim \\TTv(0,yo;w n )-y\\ 00 . T = 

where, |] . || oo;T is the uniform norm on [0, T]. If there exists a unique solution, it is 
denoted by 7Ty(0,j/o;W). 

Theorem 2.6. Let V — (Vi,...,Vd) be a collection of locally j-Lipschitz vector 
fields on R d (7 > p) such that : V and D^V are respectively globally Lipschitz 
continuous and (7 — [p])-Hdlder continuous on M. d . With initial condition yo £ M. d , 
equation dy = V(y)dW admits a unique solution 7ry(0,yo; W). 

For a proof, see P. Friz and N. Victoir j7J, Exercice 10.56. 

For P. Friz and N. Victoir, the rough integral for a collection of (7 — 1)-Lipschitz 
vector fields V — (Vi, . . . , Vd) along W is the projection of a particular full RDE's 
solution (cf. (7], Definition 10.34 for full RDEs) : dX = $(X)dW where, 

Mi = l,...,d, Va,w £ R d , $i(w,o) = (e i ,V l (w)) 

and (ei, . . . , e^) is the canonical basis of K d . 

In particular, if y : [0, T] — > A^d(M) and z : [0,T] — > R d are two continuous func- 
tions, respectively of finite p- variation and finite q- variation with 1/p + 1/q > 1, 
the Young integral of y with respect to z is denoted by y(y, z). 

Remark. We are not developing the notion of full RDE in that paper because 
it is not useful in the sequel. The reader has just to understand that a solution of 
a full RDE is a rough path. With appropriate topologies, the definition is similar 
to Definition 2.5 As mentioned above, the reader can refer to [7], Definition 10.34 
for details. 

For a proof of the following change of variable formula for geometric rough paths, 
cf. [J, Theorem 53 : 

Theorem 2.7. Let $ be a collection of ^-Lipschitz vector fields on R d (7 > p) and 
let W be a geometric p-rough path. Then, 

-1 1 



V(s,i) £ A T , $(w t )-$(w s ) 



1 



s,t 
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Now, let state some results on 1-dimensional Gaussian rough paths : 

Consider a stochastic process W defined on [0, T] and satisfying the following as- 
sumption : 

Assumption 2.8. W is a 1- dimensional centered Gaussian process with a-Holder 
continuous paths on [0,T] (a G]0, 1]). 

In the sequel, we work on the probability space (Sl,A,W) where CI = C°([0, T]; R), 
A is the cr-algebra generated by cylinder sets and P is the probability measure in- 
duced by W on (CI, A). 

Since W is an 1-dimensional process, a natural geometric 1/a-rough path W over 
it is defined by : 

(W t -W s f /a] \ 



V(«,t)€A T ,W. lt = l,W t -W s ,.. 



[1/a] 



In particular, it is matching with the enhanced Gaussian process for W provided 
by P. Friz and N. Victoir at [7], Theorem 15.33 in the multidimensional case. 

Finally, the Cameron-Martin's space of W is given by : 

H 1 w = {h£ C°([0,T];E) : 3Z G A w s.t. Vie [0,T], h t =E(W t Z)} 

with 



A w = span{W t ;t e [0,T]} . 
Let (., ■)u\ v be the map defined on Ww x 7~Lw by : 

(M)* V =E(ZZ) 

where, 

Vf G [0,T], /i t = E(W t Z) and ^ = E(W t Z) 

with Z, Z £ Aw- 

That map is a scalar product on T-L}^ and, equipped with it, is a Hilbert 
space. 

The triplet (f2,7^,P) is called an abstract Wiener space (cf. M. Ledoux |12|). 

3. Deterministic properties of the generalized mean-reverting 

equation 

In this section, we show existence and uniqueness of the solution for equation (JT|) , 
provide an explicit upper-bound for that solution and study the continuity of the 
associated Ito map. We also provide a converging approximation for equation (|T|). 

Consider a function w : [0, T] — > K satisfying the following assumption : 

Assumption 3.1. The function w is a-Holder continuous (a G]0, 1]). 

Since to is a real- valued function, a natural geometric 1/a-rough path W over it is 
defined by : 

V(s, t) G A T , W m =(l,w t - Ws ,..., {m ^f 7 "' ) • 
We also put W = / Q ] (Id[ _t] © W), which is a geometric 1/a-rough path over 

t g [o, r] .— > (t, Wt ) 
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by [7], Theorem 9.26. 



Remark. For a rigorous construction of Young pairing, the reader can refer to 
[7J, Section 9.4. 

Then, consider the rough differential equation : 

(2) dx = V(x)dW with initial condition Xq £ R, 
where V is the map defined on R + by : 

Vx £ K+, Vt, w £R, V(x).(t, w) = (a- bx)t + ax^w. 
For technical reasons, we introduce another equation : 

(3) y t =y +a(l-p) [ y^e hs ds + w t ; t £ [0, T], y > 

Jo 

where, 7 = and 

w 4 = / & s dw s with i?t = ct(1 - f3)e Hl ~ fl)t 
Jo 

for every t £ [0, T]. The integral is taken in the sense of Young. 

The map u £ [e, oo[n- u 1 belongs to C°°([e, oo[; R) and is bounded with bounded 
derivatives on [e, oo[ for every e > 0. Then, equation (pi) admits a unique solution 
in the sense of Definition |2.5| by applying Theorem |2.6| up to the time 

Tl=mi{t£[Q,T]:y t =e} ; e £}0,y ], 

by assuming that inf(0) = 00. 

Consider also the time Tq > 0, such that t\ f T o when e — > 0. 



3.1. Existence and uniqueness of the solution. As mentioned above, Section 
2 ensures that equation ([3]) has, at least locally, a unique solution denoted y. At 
Lemma 3.2 we prove it ensures that equation ^ admits also, at least locally, a 



unique solution (in the sense of D efini tion [275]) denoted x. In particular, we show 



that x — y 1+1 e~ h - . At Proposition 3.3 we prove the global existence of y by using 
the fact it never hits on [0, T]. These results together ensures the existence and 
uniqueness of x on [0, T]. 



Lemma 3.2. Consider yo > 0. Under assumptions 1 1 . l\ and 3.1 up to the time r\ 



(e £]0,yo])i if V i s the solution of with initial condition yo, then 

x:t£ [0^]^ Xt = y^ e - bt 
is the solution of ^ on [0, r^ 1 ], with initial condition xq = y^ +1 ■ 
Proof. Consider the solution y of (|3| on [0, t*], with initial condition yo > 0. 

The continuous function z — ye~ b ( 1 ~' 3 '- is bounded by M £ > on [0,t^]. 

Since 7 > 0, the map $ : u £ [0, M s ] H> u 1+1 is derivable, bounded, with bounded 
and 7-Holder continuous derivative. Moreover, by Assumption |1.1[ (7 + l)a > 1. 
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Then, by applying the change of variable formula (Theorem 2.7 1 to z and to the 
map $ between and t £ [0, t\ ] : 

rt 



i t = zj +1 + (7 + 1) / zjdz s 
Jo 

= ?7o 7+1 + f (a- bx s ) ds + a I y]e- b ^ s dw s 







Since 7 = ^(7 + 1), in the sense of Definition 2.5 x is the solution of Q on [0, r^] 
with initial condition xq = j/ 7+1 - Q 

Proposition 3.3. Under assumptions \ and \3A\ with initial condition xq > 0, 
Tg > T and then, equation admits a unique solution Try (0, Xq] w) on [0, T], 
satisfying : 

7ry(0,x ;w) = n v (0,x ; W). 

Proof. Suppose that Tq sC T, put yo = ^ ano - consider the solution y of ([3J) on 
[0, Tg] (e €=]0, j/q] ) , with initial condition y . 

On one hand, note that by definition of t^ 1 : 
Vrl ~yt = e — Vt and 

Vrl -Vt = a(l ~ /?) y yj 7 e bs ds + w T i - w t 
for every t S [0, tJ"]. Then, since t* f Tq when e — » : 

(4) y t + o(l - /3) y y^e bs ds = w t - w T i 

for every t £ [0, Tq [. 

Moreover, since u> is the Young integral of ■& £ C°°([0, T];R+) against w, and 
w is a-H61der continuous, w is also a-Holder continuous (cf. [7], Theorem 6.8). 

Together, equality ^ and the a-H61der continuity of w imply : 

-II^IU-Hoijr^o 1 - t) a ^y t + o(l - /?) / y7 7 e fcs rfs < I^IU-Hoi^^o 1 - t) a . 



On the other hand, the two terms of that sum are positive. Then, 

(5) y t «S ||w|| q _h61;t(t 1 -t) a and 

(6) 0(1 - /?) / t/7 7 e bs ds < NIU-HoLT^o 1 - 



Since t £ [0, Tq [ has been chosen arbitrarily, inequality ^ is true for every s £ [t, Tq [ 
and implies : 

Then, from inequality (pj, necessarily : 



«(i - /^Nirwy - *r ai ds < 0(1 - 0) y y s -v s d S 

||l5||a-H61;T (tq - t)" . 



Therefore, if /3 > 1 - a, t^ £ [0, T]. 
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An immediate consequence is that : 

(J [0,r E 1 ]n[0 ) T] = [0,T]. 

e6]0,yo] 

Then, (|3| admits a unique solution on [0, T] by putting : 

y = y e on [0, t}] n [0, T] 
where, y e denotes the solution of ([3]) on [0, r*] n [0, T] for every e G]0, yo]- 

In conclusion, by Lemma 3.2 equation Q admits a unique solution 7ry(0, Xq] w) 
on [0,T], matching with y 1+1 e~ b - . □ 



Remark. Note that the statement of Lemma 3.2 is true when a = 0, and up to 
the time Tq ; equation (I2| has a unique explicit solution : 



V*e [0,7b 1 ], x t = (xj 



,1-/3 



7+1 



-hi 



However, in that case, Tq ^ [0, T] is not true in general. Then, x is matching 
with the solution of equation ^ only locally. It is sufficient for the application in 
pharmacokinetic provided at Section 5. 

3.2. Upper-bound for the solution and continuity of the ltd map. Under 
assumptions 1 1 . 1 1 and 3.1 we provide an explicit upper-bound for \\ttv (0, xo', u>)lloo;T 
and show continuity results for the ltd map : 

Proposition 3.4. Under assumptions ] 1. l\ and \3. l\ for any initial condition xq > 0, 



||7fy(0, x ;w)\ 



oo;T 



1-/3 



+ o(l- (3)e bT Xo l3 T+ 



a(6V2)(l-/3)(l + T)e h ( 1 -« T || W || 00:T 



7+1 



Proof. Consider yo = x^ ^ , y the solution of ([3| with initial condition yo and 
On one hand, we consider the two following cases : 



(1) If t< 



Vt = o(l - P) 



~"<e bs ds + w T 2 



w t . 



Then, by definition of rz a 



(7) 



Vt + a(l-p) 



1 e hs ds = yo + w t — w T 2 



Therefore, since each term of the sum in the left-hand side of equality Q 
are positive from Proposition |3.3| : 



< y t < yo + \w t - w T 2 o \. 
(2) If t > Ty o ; by definition of Ty o , y t ^ yo and then, y^ 1 ^ y^ 1 '. Therefore, 
Vo < yt < yo + a(l - P)e hT ypT + \w t - w T 2 \. 
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On the other hand, by using the integration by parts formula, for every t € [0, T], 
\w t — w T 2 | = a{l -13) f e^-^dw. 

^ a(b V 2)(1 - 0)(1 + ly^-^IHU-r. 

Therefore, by putting cases 1 and 2 together ; for every t G [0, T], 

(8) < y t < yo + a{\ - P)e hT y^T + a(b V 2)(1 - + T^-^IMUt. 

That achieves the proof because, ttv(0, xq; w) — y 1+1 e~ b ' and the right hand side 
of inequality (JsJ) is not depending on t. □ 

Remark. In particular, by Proposition 3.4 ||7ry(0, xq; w)||oo;T does not explode 
when a — > or/and b — > 0. 

We remind that if it's not specified, C"- H61 ([0, T}\ R) and C°([0,T];R) are respec- 
tively equipped with ||.|| q _h61 ; t and ||.||oo ; t- 

Proposition 3.5. Under assumptions \1 . l\ and 3.1 7tv(0, .) is a continuous map 
from R* + x C a ~ H ° l ([0, T]; R) into C°([0, T]; R). 

Proof. Consider Xq,Xq > and w 1 ,w 2 : [0,T] — > R two functions satisfying As- 
sumption |3.1| 



For i = l,2, we put j/q — (^o) 1 ^ an d = ^(z/oj wnere , 

yt e [0,T], w\ = f tisdwl 
Jo 

and, with notations of equation pj, / is the map defined by : 

^(yo» = 2/o +a(l - /3) / I^(y ,w)e bs ds + w. 
Jo 

We also put : 

T 3 = mi{ S elO,T}:yl=y 2 s }. 
On one hand, we consider the two following cases : 

(1) Consider t g [0,t 3 ] and suppose that z/q ^ j/q. 

Since y 1 and y 2 are continuous on [0, T] by construction, for every s S [0, r 3 ], 
^ yl? and then, 

(//') " (//;) ' • 

Therefore, 

\Vt -Vt | = 2/t 

- Ifg - 2/o 2 + T - (rf)"> + wl - w 2 t 

Jo 

< IZ/0 ~ yol + II* 1 - W 2 ||oo;T. 

Symmetrically, one can show that this inequality is still true when y^ ^ j/q. 

(2) Consider t e [r 3 ,T], 

T 3 (t) =sup{se [r 3 ,i] :7/i=y 2 } 
and suppose that y\ ^ . 
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Since y 1 and y 2 are continuous on [0, T] by construction, for every s £ 
[r 3 0),i], y] > y 2 s and then, 

Therefore, 
\Vt ~Vt \ = Vt - Vt 

3(t) 

sc 2||iS 1 -* 2 ||o 0; t. 

Symmetrically, one can show that this inequality is still true when y\ ^ y\. 

On the other hand, by putting these cases together and since the obtained upper- 
bounds are not depending on t : 

h 1 - y 2 \\oo, T < \yl -yl\ + lT a \\w x - tD 2 || q _ H 61;t. 

Then, I is continuous from W + x C Q - H51 ([0, T}; R) into C°([0,T];R). 



For any function w : [0,T] — > M. satisfying Assumption 3.1 from Lemma 3.2 and 
Proposition |3.3| : 



Tr v (0,x ; w) 



•J7+1 



Moreover, by [7], Proposition 6.12, y(tf, .) is continuous from C a " H51 ([0, T]; R) into 
itself. 



By composition, that achieves the proof. 



□ 



Remark. Classical continuity results for the Ito map (cf. [7], Chapter 10) don't 
work here because they provide an upper-bound for \\y — ?/ 2 ||oo;T depending on : 

min y} and min y f . 

te[o,T] te[o,T] 

With similar ideas, the following proposition extends that continuity result a little 
bit. 

Consider 

C a . T = {AId [0 , T] ;AeR} x C Q - H51 ([0,T];IR), 
equipped with ||.||a,T such that : 

|M|c*,T = ||AId[ 0iT ]||(X>;T + ||w|| Q _H61;T 
= |A|T + ||w|| Q -H61;T 

for every tp = (AId [0 , T] , w) where, AeR and w E C Q - H51 ([0, T]; R). 
We also consider : 

(9) x t (X) = x + f \a-bx s (A)] Xds + a f x s (X)dw s ; t e [0,T], x > 
for any a, 6, A > 0. 



By applying Proposition 3.3 to ([9j) with a = Aa and 6 = A6, that equation ad- 
mits a unique solution tt\/(0,xq; fx,w) on [0)^1> satisfying : 

tt v (0, x ; <P\,w) = 7iv A (0, x ; iu) 
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where, ipx.w — (Aldjo^i, w) and, V\ and V are two maps respectively defined on 
R+ by : 

Vi £ R + , Vt, w e R, V x (x).(t, w) = (a- bx)Xt + ax w and 

Vx e R+, V/i, met, F(ai).(7i,, to) = (a - Sx)/i + cra'V 

Remark. Notation Ttv(0,xo;(fix,w) is justified, because it is also matching with 
ir v (0,x Q ; W x ) where, Wx = 5 [p] (AId [0 , T ] ©W). 



Existence and uniqueness of 7Ty(0, £o; Wa) is given by Proposition 3.3 via the im- 
mediate argument mentioned above. However, existence, uniqueness and other 
properties of 

tt v [0, x ; (h, w)] = n v [0, a; ; <%] (/i © W)] 
for h € C p_var ([0,T];R) such that a + 1/p > 1 will be studied in another paper. 

Now, let study the continuity of 7ry(0,xo; .) on the following subset of C a ^T '■ 
C+ T - {AId [0 , T ];A G M;} x C Q - Hal ([0,21; M) . 

Proposition 3.6. Under assumptions \ 1. l\ and \3. l\ nv{0, xq', ■) is a continuous map 
from C+ T into (7°([0, T); R). 

Proof. Let (^o = (AoId^T]: be an element of C^ T , such that = 0. 
On one hand, consider yo ~ Xq~ and the map J defined by : 

•%) =y +Aa(l-/3) / J7^)e ASs ds + w 
Jo 

for every y> = (AId[ , T ] , w) € C+ T . 

Then, for every < G [0,T], 

■Mp) - Jt(yo) = S(l - /3) / Xe xis [J-T( V ) - J7 7 (^o)] da + 

Jo 

a(l-/3) / (Ae ASs -A e A ° 6 V7 7 (w>)^ + ^-* t °- 
Jo 

With arguments of Proposition |3.5| : 

\\J(<p) - J^o)||oo;T < a(l - (3)T\\Xe xk - A e A °HUT||</~ 7 (<A))lloo;T + 

2T Q ||^-7i) |U_H61;T- 

Therefore, 

\\A<p)-J(<po)\L tT ^^o. 

¥>™>-¥>0 

In other words, J is continuous from T into C°([0, 2]; R). 



On the other hand, for any function w : [0,T] — > R satisfying Assumption 3.1 
and any A > 0, from Lemma 3.2 and Proposition |3.3| : 



tt v (0,x O ] (px, w ) = nv x (0,x o ;w) 

= e- x ~»-,r +1 {\id [0 , Th y($x,w)} 

where, d x = <r(l - ( S)e AS ( 1 "«' and ^ = (AId [0 , T ], 

Moreover, by [7], Proposition 6.12, y is continuous from 

C<i-H61q 0)T ] ;R ) x c* Q - H51 ([0,T];R) into C"- H61 ([0, T]; 
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By composition, that achieves the proof. 



□ 



Remark. As mentioned above, by Proposition 3.4 ||7ry [0, xq; (AIcWxi) ^)]||oo;T 
does not explode when A — > 0. Precisely, there exists at least one sequence (A n , n € 
N) of elements of and one continuous function L : [0, T] — > E + such that : 



lim A„ = and lim \\ttv \0, Xq', (A„Idr t] > w )] ~ L\ 



oo;T 



0. 



3.3. A converging approximation. In order to provide a converging approxima- 
tion for equation Q , we first prove the convergence of the implicit Euler approxi- 
mation (y n , n G N* ) for equation ([3| : 



(10) 



Va 



yo > o 



Vk+i =Vk + 



a(l-/3)T 



where, for n e N* , = kT/n and k < n while y£ +1 > 0. 



The following proposition shows that the implicit step-n Euler approximation y n is 
defined on {1, . . . , n} : 

Proposition 3.7. Under assumptions \ 1. l\ and 3.1 equation (10) admits a unique 
solution (y n ,n € N*). Moreover, 

ViiGM*, Vfc = 0,...,n, yl >0. 

Proof. Let / be the function defined on x E x by : 

Va £ K, Vcc, b > 0, / (x, a, b) = x — bx^ 1 — a. 

On one hand, for every a £ K and b > 0, f(.,a,b) £ C°°(R+;K) and for every 
x > 0, 

d x f {x, a, b) = 1 + bjx'^ +1) > 0. 
Then, f(.,a,b) increase on M.* + . Moreover, 

lim / (a;, a, 6) = — oo and lim / (x, a, b) = oo. 

Therefore, since / is continuous on x K x Ml : 

(11) Va € M, V6 > 0, 3\x > : f(x,a,b) = 0. 



On the other hand, for every nfiN*, equation (|10| can be rewritten as follow 
(12) 



a(l — 0)T ht n „ 



n 



= 0. 



In conclusion, under Assumption |l.l| by applying ( |11[ ) to equation (12) for each 
integer k between and n, y n is defined on {1, . . . , n}. 



Necessarily, y^ > for k = 0, . . . , n. 
That achieves the proof. 

For every n £ N* , consider the function y n : [0,T] 



a 



k=0 



Vk 



Vk+l Vk /, _ , n \ 
t fc+l *; 



such that 

1 [tjf,tff +1 [(*) 



for every t€ [0,T]. 



The following lemma provides an explicit upper-bound for (n, () e N* x [0, T] i-t y™. 
It is crucial in order to prove probabilistic convergence results at Section 4. 
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Lemma 3.8. Under assumptions \l . l\ and \3. i] : 

sup HinL-r < yo + a(l-(3)e bT y^T + 

nGN* 

a(6V2)(l-^)(l + T)e fc ( 1 -« T || W || 00;T . 
Proof. Similar to the proof of Proposition |3.4| 



First of all, by applying ( 10 I recursively between integers I < k n and a 
change of variable : 

(13) Vl - VT = a{l ~ P)T E (»JV e bt " + w t n - w t? . 

i=l+l 

Consider n e N* and 

k yo = max {A; = 0, . . . , n : y k < y } ■ 
For each k = 1, . . . , n, we consider the two following cases : 
(1) If k < k yo , from equality (131 : 

k 

a{l-p)T ^ ( „ A _ 7 bt? 

Uk vo yk - 

Then 

(14) vi + £ forr 7 e " tr = < + 0* - ffi. 



i=k+l 



o(l-/3)T ^ 



i=fc+l 



Therefore, since each term of the sum in the left-hand side of equality ( 14) 
are positive from Proposition |3.7| : 



k 

n 

i=k+l 

< yo + \m? - w t n | 

because y? < y - 

(2) If fc > fc ao ; by definition of k yo , for i — + 1, . . . , k, yf > yo and then, 
(j/™) -7 ^ y^ 1 . Therefore, from equality (13) : 

i=k yQ -\-l 

sC y + o(l - /?)e bT % 7 T + |u) t « - ^ |. 
As at Proposition |3.4| : 



sup y 7 t l ^ max y£ 

t£[0,T] k=0,...,n 

(15) < 2/o + o(l - ^)e fcT y - 7 T + <r(6 V 2)(1 - /3)(1 + iV^^IM 

That achieves the proof because the right hand side of inequality (15) is not de- 
pending on n. □ 

With ideas of A. Lejay [13], Proposition 5, we show that (y n ,n € N*) converges 
and provide a rate of convergence : 

Theorem 3.9. Under assumptions \l . l\ and \3. l\ (y n ,n £ N*) is uniformly converg- 
ing to y, the solution of equation ftSw for initial condition y , with rate n" 01 " 1111 ^ 1 ' 1 '. 
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Proof. It follows the same pattern that Proof of [13], Proposition 5. 

Consider n G N*, t G [0, T] and y the solution of equation ^ with initial con- 
dition y > 0. Since (t2;k = 0, ...,n) is a subdivision of [0, T], there exists an 
integer k ^ n - 1 such that t G 

First of all, note that : 

(16) \y? -yt\< \vt -vl\ + \Vk ~ 4\ + VI - vA 

where, zf = yt™ for i — 0, . . . , n. Since y is the solution of equation (|3j), z£ and 
satisfy : 

4+i = 4 + fl(1 ~ /3)r (zlVi)' 7 ^ 1 + tf* s+1 - ms + 4 

where, 

4 = o(i - /3) r +1 (y s " 7 e bs - 

In order to conclude, we have to show that \y% — zjE? | is bounded by a quantity not 
depending on k and converging to when n goes to infinity : 

On one hand, for every [u, v) G A^, 



| e ^-7_ e f-y-7| = 



e bv yZ - e hu yl 



7 7 



^ Tz^(z bv \yl-yv\ + \yvV\e bu -e b v\ 

WuVv I 



V-'WLt (\\y\C&\v - u\°™M + b\\y\\l ;T \v - u\) 



because s G M+ M> s 7 is 7-Holder continuous with constant 1 if 7 G]0, 1] and 
locally Lipschitz continuous otherwise, y admits a strictly positive minimum and 
is a-H61der continuous, and s G [0, T] 1— > e bs is Lipschitz continuous with constant 
be bT . In particular, if |u — u\ S% 1, 

|e b V 7 " e b "y- 7 | < e bT \\y^\\l lT (|M|£. Hel;T + %||L ; t 

where /1 = min(l,7). 

Then, for i = 0, . . . , k, 

< a(l - /3) lyjV 8 - y^/^ds 



^a(l-^)||eS- 7 |L u . HH1: T / ' (t^-s^ds 



a/i-Hol;T y n 
a 0- ~ P) \\ e b. „-7|| __>_ 



(171 < h 

^ ^ aM + l 'I ^ MaM-HaiiTjjQMi+l- 

On the other hand, for each integer i between and k — 1, we consider the two 
following cases (which are almost symmetric) : 

(1) Suppose that yf +1 > z£ +1 . Then, 
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Therefore, 

llJi+l ~ z i+l\ — Ui+1 ~ z i+l 

= yf - z? + a{1 -J )T e bt ^ [{y? +1 )-< - (zJVi)" 7 ] 

< |y»-Z«| + | £ ™|. 

(2) Suppose that zf +1 > yf +l . Then, 

fei)" 7 -fei)" 7 <o. 

Therefore, 

\ z i+l ~ Di+l\ = z i+l ~ Ui+1 

= Z? _ y n + a 0^M e K +1 [( ^ +i) - 7 _ (y n +i) - 7] + £ n 

< \yV; - Z?\ + W\. 
By putting these cases together : 

(18) Vi = 0, . . . , k - 1, k" + i - tf+il < I*" - 2/"l + 



By applying (18) recursively from k — 1 down to : 

k — 1 

< |2/o-zo|+X>?l 

i=0 

(19) ^ 0(1-^) ^+111 6 - 7 || J >Q 



because yo = zq and by inequality (17). 



Moreover, from inequality (191, there exists N 6 N* such that for every integer 
?i > AT, 

K+i " 4+1 1 < . max \y? - *?| < m y 

2— l,...,n 

where, 

1 

m y = - mm j/ s . 

2 s£[0,T] 

In particular, 

2/fe+i > z k+i ~ m y^ my. 

Then (y£ +1 )~ 7 < ra~ 7 , and 



i j.n 4-n 

L k+1 L k 

^ [o(l - /3)Te &T m" 7 + |H| a _ H 6i;T] ► °- 



In conclusion, from inequality (16 1 : 

(20) |y? - y t |< [o(l ~ /3)Te bT m-; 7 + |HU_ H 6i ; T + |MU-Ho1 ; t] ^ 

aj 1 ~ P) T a»+1 lib - — y II _i 

a/J + 1 " lla/<-H61;T n op n^co 



That achieves the proof because the right hand side of inequality ( 20 ) is not de- 
pending on t. □ 
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Finally, for every n g N* and t £ [0,T], consider cc" = e h *(y™) 7+:L . 

The following corollary shows that [x n ,n £ N*) is a converging approximation 
for x = tt(0,Xq;w) with x > 0. Moreover, as the Euler approximation, it is just 
necessary to know xq, w and, parameters a, 6, a > and /3 > 1 — a to approximate 
the whole path x by x n : 



Corollary 3.10. Under assumptions and 3.1 (x n ,n £ N*) is uniformly con- 
verging to x with rate 7i _Qmin ( 1 >')')_ 

Proof. For a given xq > 0, we shown that x = e~ b y 7+1 is the solution of equation 
^ by putting yo = Xq~^ , where y is the solution of equation ^ with initial con- 
dition ?/o- 



From Theorem 13.91 : 

Wx-x^t ^ C\\y - y n \\oo;T 

sC C [a(l - l3)Te bT m-~> + \\w\\ a . m i.T + IMU-HoLt] ^ 

C ail ~^ T^ +1 \\e b -y-n >0 

a/X + 1 11 lla/i-Hol;T n aji „^oo 

where, C is the Lipschitz constant of s n- s 7+1 on 

0, \\y\\oo-T + SUp ||y"||oo;T 
ngN* 

Then, (x n ,n g N*) is uniformly converging to x with rate n" Qmm ( 1,7 ^. □ 

Remark. When a > 1/2 ; (3 > I- a > 1/2 and then 7 > 1. Therefore, (i",nef) 
is uniformly converging with rate n~ a < n 1-2 ". In other words, the approximation 
of Corollary |3. 10| converges faster than the classic Euler approximation for equations 
satisfying assumptions of |13) . Propositions 5. 

4. Probabilistic properties of the generalized mean-reverting 

equation 

Consider the Gaussian process W and the probability space (n,,A,V) introduced 
at Section 2. Under Assumption |2.8| almost every paths of W are satisfying As- 
sumption |3T| Then, under assumptions |1 . 1| and 2.8 results of Section 3 hold true 



for Try (0, xq; W), with deterministic initial condition xq > 0. 

This section is essentially devoted to complete them on probabilistic side. In par- 
ticular, we prove that 7Tv{0, Xq; W) belongs to L p (f2) for every p ^ 1. We also show 
that the approximation introduced at Section 3 for ity (0, xq; W) is converging in 
L p (n) for every p ^ 1. 

Remark. Since W is an 1-dimensional process, as mentioned at Section 2, there 
exists an explicit geometric 1/a-rough path W over it. That explains why Assump- 
tion 



2.8 is sufficient to extend deterministic results of Section 3 to ttv(0,xq; W). 



4.1. Basic probabilistic properties. In this subsection, we show what happens 
to X = 7Ty(0,a;o; W) when W has stationary increments or satisfies a self-similar 
property : 

Proposition 4.1. Assume that W satisfies Assumption \2.8\ and there exists h > 
such that : 

w, +h - w h = w. 
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Under Assumption \l.l\ for every deterministic initial condition xq > 0, 
TT V]O , t+h {0, x ; W) = 7t v . 0>t (0, X h ; W) 

for every t £ [0, T] . 

Proof. By Proposition |3.3| X has almost surely continuous and strictly positive 



paths on [0, T]. Then, by Theorem 2.7 applied to almost every paths of X and to 
the map u i-> m 1 ~' 9 between and t £ [0, T] : 

X\- & = x\-^ + {1-0)1 X-?(a - bX u )du + a(l - f3)W t . 
Jo 

Therefore, X X 7^ = Z(h) where, 

Z t (h) = X 1 -? + (1-/3) f Z-\h) [a - bZZ+\h)] du + (7(1 - p)W t ; * G [0,T] 
Jo 

because W.+h — Wh = W. 



In conclusion, by applying Theorem 2.7 to almost every paths of Z(h) and to 
the map u i-> w 7+1 : 

X t+h -X h = (a- bX u+h ) du + a X? +h dW u 
Jo Jo 

for every t G [0,T]. □ 



Proposition 4.2. Assume that W satisfies Assumption 2.8 and there exists h > 
such that : 

Ve > 0, W e . = e h W. 
Under Assumption ]! . 1\ for every deterministic initial condition xq > 0, 

7TV;0,et (0, X \ W) = ^V e>h ;0,t (0, X ] W) 

for every t G [0, T] and e > 0, with : 

Viz; G R+, Vt,w G K, V £< h(x).(t,w) = e(o - + crs h x p w. 
Proof. By Proposition |3.3| X has almost surely continuous and strictly positive 



paths on [0, T]. Then, by Theorem 2.7 applied to almost every paths of X and to 
the map u i-> w 1- ' 9 between and t G [0, T] : 



x\- & = x\- & + (1-/3) X-f>(a - bX u )du + a{l - f3)W t . 
Jo 

Therefore, for every e > 0, X\~$ = Z(e) where, 

Z t {E)=x\- p +e(l-p) f Z-~I(e) [a-bZZ +1 (e)]du + e h a{l~ $)W t ; t € [0,T] 



because W e = e h W. 



In conclusion, by applying Theorem 2.7 to almost every paths of Z(s) and to the 
map u I—)- u 7+1 : 



X e t = Xq 



el (a- bX £U ) du + ae h f X? u dW u 
Jo Jo 

for every t G [0, T] and e > 0. □ 



Remark. Typically, mean-reverting equations driven by a fractional Brownian 
motion are concerned by propositions |4.1| and |4.2| 
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4.2. Integrability and convergence results. Consider the implicit Euler ap- 
proximation (Y n , n e N* ) for the following SDE : 

r-t 

Yt = Vo + a{ 



where, 



i(l - /3) f Y-~>e bs ds + W t ; t G [0,T], y > 
Jo 

W t = I $ s dW s and d t = <r(l - /3)e b(1 -' 9)t 
Jo 

for every i € [0,T]. 

Proposition 4.3. Under assumptions \l.l\ and \2.£\ for every deterministic initial 
condition x > 0, 

(1) ||7Ty(0, Xq; W)\\oo-,t belongs to L p (fl) for every p ^ 1. 

(2) For every p ^ 1, 

™p ||X"IL ;T G 

where, for every n G N*, -X" n = e — b '(l r ™) T+1 2/0 = ^0 ■ 
Proof. On one hand, by Proposition |3.4| and Fcrnique's theorem : 



||7iv(o,flso;W)|| 



for every p ^ 1. 



+ a(l - /3)e 6T z^T+ 
< 7(6V2)(l- i 8)(l + r)e 6 ( 1 -« r ||W|| 00 .r 



7+1 



On the other hand, by Lemma |3.8| and Fernique's theorem : 
sup imioo-T < yo + a(l-(3)e bT y^T + 

a(b V2)(l - /3)(1 + lOe^^^HWlU-r G 
for every g ^ 1. Then, by putting q = (7 + l)p for every p ^ 1, 

sup phi T eL*(n). 

neN* 



□ 



Corollary 4.4. Under assumptions 1.1 and 2.8, for every deterministic initial 
condition xq > 0, (X n ,n G N*) is uniformly converging to 7ry( u i Xo', w) in L p (fl) 
for every p ^ 1 . 



Proof. By Corollary pUP] 



p^-frvKO.ZojWOIL^ 



P-a.s. 



> 0. 



Then, by Proposition 4.3 and Vitali's convergence theorem, (X n ,n G N*) is uni- 
formly converging to 7ry(0, xq; W) in L p (f2) for every p ^ 1. □ 



Remark. Note that Proposition 4.3 is crucial to ensure this convergence in L p (fl) 



for every p ^ 1. Indeed, inequality (20 1 doesn't allow to conclude because it is not 
sure that \\e b -Y~-<\\ a ^6Y,T G i 1 ^)- 

For example, let simulate (X n ,n G N*) with T = 1, xq = 1, a = 6, 6 = 4, 
a = 1, P = 0.8 and a fractional Brownian motion B H with Hurst parameter 
i? G {0.3,0.8} : 
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4.3. A large deviation principle for the generalized M-R equation. We 

establish a large deviation principle for the generalized mean-reverting equation (as 
P. Friz and N. Victoir at (7], Section 19.4). 

First of all, let's remind basics on large deviations (for details, the reader can 
refer to 0). 

Throughout this subsection, assume that inf(0) = oo. 

Definition 4.5. Let E be a topological space and let I : E [0, oo] be a good 
rate function (i.e. a lower semicontinuous map such that {x € E : I(x) A} is a 
compact subset of E for every A ^ 0). 

A family (fi £ ,e > 0) of probability measures on (E,B(E)) satisfies a large devi- 
ation principle with good rate function I if and only if for every A 6 B(E), 

-I(A°) sc lim^log \ii e (A)] < fim^oelog [/x £ (A)] < -1(A) 

where, 

VAeB(E), 1(A) = inf I(x). 

Proposition 4.6. Consider E and F two Haussdorff topological spaces, a continu- 
ous map f : E F and a family (fJ, e ,e > 0) of probability measures on (E,B(E)). 

If (/i £ ,e > 0) satisfies a large deviation principle with good rate function I : E — > 
[0, oo], then (/i £ o / , e > 0) satisfies a large deviation principle on (F, B(F)) with 
good rate function J : F — > [0, oo] such that : 

J(y) = inf {I(x);x € E and f(x) = y} 

for every y G F. 



20 



NICOLAS MARIE 



That result is called contraction principle. The reader can refer to [2], Lemma 4.1.6 
for a proof. 

Consider the space C°' Q ([0, T]; K) of functions ip G C Q " H51 ([0, T]; E) such that : 

hm uJS) = o with w „(<n = sup 

\t - s\ ^ j 

for every 8 > 0. 

In the sequel, C 0,a ([0, T]; M.) is equipped with ||.|| q -H61;T an d the Borel cr-field gener- 
ated by open sets of the a-H61der topology. The same way C° ([0, T] ; M) is equipped 
with ||.|joo ; T and the Borel cr-field generated by open sets of the uniform topology. 

Now, suppose that W satisfies : 
Assumption 4.7. There exists h > such that : 

Ve > 0, W s . = e h W. 

Moreover, C C°- a ([0,T];R) and (C°' a ([0,T];lR),K^,P) is an abstract Wiener 
space. 

Remarks : 

(1) The notion of abstract Wiener space is defined and detailed in M. Ledoux [12J 

(2) Typically, the fractional Brownian motion with Hurst parameter H > 1/4 
satisfies Assumption 4.7 (cf. [17], Proposition 4.1). 

Consider the stochastic differential equation : 

(21) X t = x + ±J (a- bX s ) ds + ^ X^dW a ; t G [0, T] 

where, xq > is a deterministic initial condition, a, b, a, 8 > and j3 G]0, 1] satisfies 
Assumption |1.1| 

Under assumptions |1.1| and 2.8 by propositions 3.3 and 4.3 equation (21 1 ad- 
mits a unique solution belonging to L p (fl) for every p ^ 1. 

Moreover, under Assumption |4.7| by Proposition |4.2| : 

(22) X st =x + £ -J^ (a - bX ss ) ds+-^ J X? s dW s 
for every t G [0, T] and e > 0. 

In the sequel, assume that <5 = e. Then, X s , satisfies : 

X e .=n v (0,x ;eW) 
where, V is the map defined on K + by : 

G K+, W, Kiel, V(a;).(t, w) = (a - 6x)i + ctx'V 
Let show that (X e ., e > 0) satisfies a large deviation principle : 



Proposition 4.8. Consider xo > 0. Under assumptions \l . l\\2.8\ and Jf..1 (X e .,e > 
0) satisfies a large deviation principle on C°([0,T];K) with good rate function J : 
C°([0,T];R) [0,oo] defined by : 

Vy G C°([0,T];]R), J(y) = inf {l(w); w G C Q '°([0, T]; M) and y = 7f v (0,a:o; wj) 
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where, 

1(W) -\ ocifwfU^ 
for every w € C°' Q ([0, T]; M). 



3.5 



Proo/. Since C°' a ([0, T];R) C C Q - H51 ([0,T];R) by construction, Proposition 
implies that 7ry(0, xo; .) is continuous from 

C°' a ([0,T];M) into C°([0,T];R). 

On the other hand, under Assumption |4.7| by M. Ledoux [T^], Theorem 4.5 ; 
(eW,e > 0) satisfies a large deviation principle on C°' Q ([0, T]; R) with good rate 
function /. 

Therefore, since X e _ = 7ry(0, a^o; eW) for every e > 0, by the contraction principle 
(Proposition 4.6 1, (X e ,e > 0) satisfies a large deviation principle on C°([0, T];R) 
with good rate function J. □ 

5. A GENERALIZED MEAN-REVERTING PHARMACOKINETIC MODEL 

We study a pharmacokinetic model based on a particular generalized mean-reverting 
equation (inspired by K. Kalogeropoulos et al. 

In order to study the absorption/elimination processes of a given drug, the fol- 
lowing deterministic mono-compartment model is classically used : 

(23) C t = £ (^^e- K " 8 - K e c)j ds:te [0, T] 

where : 

• A > is the dose administered to the patient at initial time. 

• v > is the volume of the elimination compartment E (extra-vascular 
tissues) . 

• K a is the rate of absorption in compartment A. If the drug is ad- 
ministered by rapid injection, an IV bolus injection, it is natural to take 
K a = 0. 

• K e > is the rate of elimination in compartment E, describing removal of 
the drug by all elimination processes including excretion and metabolism. 

• Ct is the concentration of the drug in compartment E at time t € [0, T]. 
Remark. About deterministic pharmacokinetic models, the reader can refer to Y. 
Jacomet |5]. 

Recently, in order to modelize perturbations during the elimination processes, sto- 
chastic generalizations of ( 23 1 has been studied : 

Ct = L {^r e KaS ~ KeCs ) ds + l a (s ' Cs) dB *' te t°' T i 

where, B is a standard Brownian motion and the stochastic integral is taken in the 
sense of Ito. For example, in K. Kalogeropoulos et al. |TT] : 

C t = j (^^e- K « s - K e C s ^j ds + aj^ C?dB s ; t € [0, T] 

with a > and (3 £ [0, 1]. 

However, these models aren't realistic (cf. M. Delattre and M. Lavielle [3]), be- 
cause the obtained process C is too rough. 
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Since probabilistic properties of Ito's integral aren't particularly interesting in that 
situation, if the drug is administered by rapid injection, C could be the solution of 
equation with Co — Aq/v, a — and b = K e . 



In order to bypass the difficulty of the standard Brownian motion's paths rough- 



ness, one can take a Gaussian process W satisfying Assumption 2.8 with a close 
to 1. Typically, a fractional Brownian motion B H with a high Hurst parameter H 
(cf. simulations below). 



Precisely 
(24) 



C t = — - K e [ C s ds + a [ CgdWs 
v Jo Jo 



where the stochastic integral is taken pathwise, in the sense of Young. Moreover, 
since a = 0, we shown at Section 3 that until it hits zero, the solution of equation 



(24 1 is matching with the process X defined by 

7+1 



Vt e 



Ao 

v 



1-/3 



w t 



- K ' 1 with W t = <t(1 - P) / e^-ft'dWt 



It is natural to assume that when the concentration hits 0, the elimination process 
stops. Then, we put C — X\^ T i^ T y where T > is a deterministic fixed time. 

For example, let simulate that model with Ao — v, K e — 4, a — 1, (3 = 0.8 
and a fractional Brownian motion B H with Hurst parameter H € {0.6, 0.9} : 



M- R model v.s. deterministic model, H = 0.9 



M- R model v.s. deterministic model, H = 0.6 
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Figure 2. GM-R model v.s. deterministic model currently used 



Even with H = 0.6, the solution seems not to be realistic. With H — 0.9, there is 
just little perturbations in the elimination process. 

In the sequel, we also consider the process Z = X 1 ^ 13 . Its covariance function 
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2:5 



is denoted by cz- 



For clinical applications, parameters K e , a and j3 have to be estimated. Con- 
sider a dissection (t , . . . ,t n ) of [0, T] for n G N*. We also put Xi = X t . and 
z.j = Z ti for i = 0, . . . , n. The following proposition provides the likelihood func- 
tion of (x\, ...,x n ) which can be approximatively maximized with respect to the 
parameter 9 = (K e , a, /3) by various numerical methods (not studied in this paper) : 



Proposition 5.1. Under assumptions \l . 1\ and 2.8 the likelihood function of [x\, . . . , x n ) 
is given by : 



L(6;xi,...,x n ) 



2"(1 — P) U ~i-x 1 >0,...,x n >0 TT -p x 

(2tt)»/2 ^det [T{K ei a,P)\ * l 



exp 



■~{T-\K e ,a,p) 



-0 _ C l-P e -K e (l-t 



(x\- p -Cl- p e- K ^-^ 



whe 



a 2 (6) = Var(z u ...,z n ) and T{9) 



a\{9) 



Cz(tn,h) 



cz(h,t n ) 



Proof. Since W is a centered Gaussian process as a Wiener integral against W ; 
Z\, . . . ,z n are n Gaussian random variables. We denote by fi n {6; .) the natural 
density of (z\, ... , z n ) with respect to the Lebesgue measure on (R n , B(R")). 

Consider an arbitrary Borel bounded map if : K" — > K. By the transfer theo- 
rem : 

E{cp(x u ...,x n )} =E[<p(\z 1 \T+ 1 ,...,\zZ+ 1 \)] 

= 2" / <p(a} +1 , . . . ,a1 +1 )fi,...,n{&;a 1 , . . . ,a n )da\ . . . da n 

by reduction to canonical form of quadratic forms. 

Put Ui = aj +1 for dj € and i = 1, . . . , n. Then, 



(ai, . . . ,a„) = (uf 



,in +1 ) and |J(ui, . . . ,u n ) 



(7- 



where, J(ui, . . . , u„) denotes the Jacobian of ; 



( Ul ,..., Un ) g (»;) n 1— > 

By applying that change of variable : 

2™ 



E[ip(x 1 ,...,x n )} = 



dui . . . du n tp(ui, . . . , u n ) x 



(7 + I) 1 

/i,...,n(fl;«?* r ,...,«^ r )n i 
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Therefore, ¥t Xl x \(6; dui, . . . , du n ) = L{6; u\, . . . , u n )dui . . . du n with 



2" J— J— TT 1 

L(6;ui,...,u n ) = ^ - ^ n fi,...,„(0; uj +1 uZ +1 ) [[ u % 1+1 l Ul 



>0....,u„>0 



2"(1 — /3)"1 Mi >q,..., m „>q -pr -f) 

(27r)»/Vdet [T(K e ,a,f3)]l\ Ul 

/u\~ p - Cl- p e- K ^-P^ 



exp 



\ u l-f> _ C l-P e -K e (l-P)t„ 



□ 



Finally, consider a random time r £ [0, Tq A T] and a deterministic function F : 
R + —> M. satisfying the following assumption : 

Assumption 5.2. The function F belongs to C 1 (M.+ ;M.) and there exists (K, N) £ 
R* + x N* such that : 

Vr £ M +; |F(r)| K(l + r) N and \F(r)\ sC K(l + r) N . 

Let show the existence and compute the sensitivity of f T (x) = E[F(C*)] to varia- 
tions of the initial concentration x > in compartment E. 



Proposition 5.3. Under assumptions 1.1\ 2.8 and \5.2\ the function f T is derivable 



on R!j_ and, 

\/x > 0, fr(x) = x- E \e- K " T F{C x T ){x 1 - 13 + W T )^ . 
Proof. First of all, the function x £ M> is almost surely C 1 on and, 

Vz > 0, d x C x T = x-V (x x ~^ + e~ K * r . 

Consider x > and s £]0, 1]. 

On one hand, since F belongs to C 1 (IR+;IR), from Taylor's formula : 



F(C* +£ ) - F(C*) 



F{C* +ee )d x C* +ee d6 



< sup K{\ + \\C BJ ^ e \\ ooi r) N \d m C^ 4>e \ 
0e[o,i] 



by Assumption |5.2| 

On the other hand, since 8, e £ [0, 1] 



(25) 
and 
(26) 



\\C x+ee \U; T < Ux + l) 1 -^ + \\W\\ 00 , 7 



\d x C* +ee \ sc x-P \{x + I) 1 - + \\W\l 



7+1 



By Fcrnique's theorem, the right hand sides of inequalities (25) and (26) belong to 



L P (Q) for every p > 0. Moreover, these upper-bounds are not depending on 9 and e. 
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Therefore, by Lebesgue's theorem, f T is derivable at point x and, 
f T (x) = :ir% \er K " T F(C^)(x 1 - f3 + WrV ■ 



a 



There is probably many ways to use that result in medical treatments. For ex- 
ample, assume that f T (x) modelize a part of patient's therapeutic response to the 
administered drug. Proposition |5.3| provides a way to minimize the initial dose for 
an optimal response. 

Remarks : 

(1) By the strong law of large numbers, there exists an almost surely converging 
estimator for that sensitivity. 

(2) For any x > 0, one can show the existence of a stochastic process h x defined 
on [0,T] such that f T (x) — E[F(C x )5(h x )] where, S denotes the divergence 
operator associated to the Gaussian process W . Then, F has not to be 
derivable anymore by assuming that F G L 2 (H$+). It is particularly useful 
if F is not continuous at some points. 



We don't develop it in that paper because the Malliavin calculus frame- 
work has to be introduced before. To understand that idea, please refer to 
E. Fournie et al. |B| in Brownian motion's case and N.M. |10| . 
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